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Abstract 

Dynamic allocation of resources to the best link in large multiuser networks offers considerable 
improvement in spectral efficiency. This gain, often referred to as multiuser diversity gain, can 
be cast as double-logarithmic growth of the network throughput with the number of users. In 
this paper we consider large cognitive networks granted concurrent spectrum access with license- 
holding users. The primary network affords to share its under-utilized spectrum bands with the 
secondary users. We assess the optimal multiuser diversity gain in the cognitive networks by 
quantifying how the sum-rate throughput of the network scales with the number of secondary 
users. For this purpose we look at the optimal pairing of spectrum bands and secondary users, 
which is supervised by a central entity fully aware of the instantaneous channel conditions, 
and show that the throughput of the cognitive network scales double-logarithmically with the 
number of secondary users (N) and linearly with the number of available spectrum bands (M), 
i.e., M log log iV. We then propose a distributed spectrum allocation scheme, which does not 
necessitate a central controller or any information exchange between different secondary users 
and still obeys the optimal throughput scaling law. This scheme requires that some secondary 
transmitter-receiver pairs exchange log M information bits among themselves. We also show that 
the aggregate amount of information exchange between secondary transmitter-receiver pairs is 
asymptotically equal to A/dog A/. Finally, we show that our distributed scheme guarantees 
fairness among the secondary users, meaning that they are equally likely to get access to an 
available spectrum band. 

kewords: Cognitive radio, distributed, fairness, multiuser diversity, spectrum allocation. 



1 Introduction 

Dense multiuser networks offer significant spectral efficiency improvement by dynamically identify- 
ing and allocating the communication resources to the best link. The improvements thus attained 
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are often referred to as multiuser diversity gain and rest on the basis of opportunistically allocating 
all the resources to the most reliable link. The performance of such resource allocation scheme 
relies on the peak, rather than average, channel conditions and improves as the number of users 
increases, as it becomes more likely to have a user with an instantaneously strong link. 

The notion of opportunistic communication and multiuser diversity was first introduced [1] for 
uplink transmissions, and further developed in [2H3] for downlink transmissions. The analysis of 
multiuser diversity gain in downlink multiple-input multiple-output (MIMO) channels is provided 
in [5j[6]. In all these transmission schemes, the sum-rate capacity exhibits a double- logarithmic 
growth with the number of users. 

The recent advances in secondary spectrum leasing [7] and cognitive networks [8] suggest ac- 
commodating unlicensed users (secondary users or cognitive radios), within the license- holding 
networks and allowing them to access under-utilized spectrum bands. Among different spectrum 
sharing schemes, underlaid spectrum access [9] has received significant attention. This scheme 
allows for simultaneous spectrum access by the primary and secondary users, provided that the 
power of secondary users is controlled such that they impose limited interference to the primary 
users. 

In this paper we consider opportunistic underlaid spectrum access by secondary users and 
assess the multiuser diversity gain by analyzing the sum-rate throughput scaling of the cognitive 
network. Such analysis for cognitive networks differs from those of primary networks studied 
in [THSIEIE] in two directions. First, the transmissions in the cognitive network are contaminated 
by the interferences induced by the primary users. The existence of such interference does not 
make opportunistic communication possible by merely finding the strongest secondary link, and 
necessitates accounting for the effect of interference as well. Secondly, and more importantly, the 
uplink/downlink transmissions in the networks referenced above, require feedback from the users to 
the base station and it is the base station that dynamically decides which user(s) should receive the 
resources. Cognitive networks, in contrary, are often assumed to lack any infrastructure or central 
entity and spectrum allocation should be carried out in a distributed way. 

To address these two issues, we first focus on only examining the for the effects of interference and 
assume that the cognitive network has a central decision-making entity, fully aware of all cognitive 
users' instantaneous channel realizations. This result, providing the optimal scaling factor, presents 
an upper bound on the throughput yielded by any distributed spectrum allocation scheme. In the 
next step, we offer a distributed algorithm where the secondary users decide about accessing a 
channel merely based on their own perception of instantaneous network conditions. 

Our analyses reveal that, interestingly, in both centralized and distributed setups, the sum-rate 
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throughput scales double-logarithmically with the number of users, which is the optimal growth and 
is the same as that of centralized primary networks. Therefore, the interference from the primary 
network incurs no loss on the multiuser diversity gain of the cognitive network. 

We also examine how fairness is maintained in our distributed scheme. Generally, in opportunis- 
tic communication schemes there exists a conflict between fairness and multiuser diversity gain, as 
the network tends to reserve the resources for the most reliable links, which leads the network to 
be dominated by the users with strong links. We show that, however, in our distributed scheme, 
we can ensure fairness among the secondary users by providing them with the same opportunity 
for accessing an available spectrum band. 

The remaining part of the paper is organized as follows. In Section [2] we provide the system 
model as well as the statement of the problem. Sections [3] and H] discuss the sum-rate through- 
put scaling laws in centralized and distributed cognitive networks, respectively. Our distributed 
algorithm requires some information exchange between each cognitive transmitter and its desig- 
nated receiver. The amount of such information is quantified in Section 15. 1L As we are considering 
opportunistic type of spectrum access, it is crucial to also look at the fairness among the users. 
The discussion on the fairness is given in Section 15.21 Some remarks on the implementation of the 
distributed spectrum access algorithm are provided in Section [6] and Section [7] concludes the paper. 
In order to enhance the flow of the material, we have confined most of the proofs in the appendices. 

2 System Descriptions 
2.1 System Model 

We consider a decentralized cognitive network comprising of N secondary transmitter-receiver pairs 
coexisting with the primary transmitters via underlaid [9] spectrum access. Therefore, the primary 
and secondary users can coexist simultaneously on the same spectrum band. The primary network 
affords to accommodate 1 < M -C N secondary users and allows them to access the non-overlapping 
spectrum bands B±, ... , Bm such that each band is allocated to exactly one secondary transmitter- 
receiver pair. We assume that the secondary transmitters and receivers are paired up a priori such 
that each secondary transmitter knows its designated receiver and vice versa. We also assume 
that each secondary transmitter and receiver is potentially capable of operating on each of the M 
spectrum bands, a feature facilitated by having appropriate reconfigurable hardware. 

We assume quasi-static flat fading channels and denote the channel from the j th primary trans- 
mitter to the i th secondary receiver in the m th spectrum band (B m ) by /i™ £ C and denote the 
channel between the i th secondary transmitter-receiver pair in the m th spectrum band (B m ) by 
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g™ 6 C. Let x?(t) and xf(t) represent the transmitted signals by the i th primary transmitter and 
the i th secondary transmitter, respectively. We assume that there might be a group of active pri- 
mary users on each spectrum band B m and define the set B m such that it contains the indices of 
such users. If the n th secondary pair transmits on B m , then the received signal at the n th secondary 
receiver is given by 

Vn = V«< + £ X~'C/j + C (1) 

jeB m 

where z™ E CJ\f(0, Nq) is the additive white Gaussian noise at the n th receiver. In a non- 
homogeneous network, the users experience different path-loss and shadowing effects, which we 
account for by incorporating the terms {jij} and {rji}. Also, we assume that the primary and 
secondary transmitters satisfy average power constraints P p and P s , respectively, i.e., E[|x^| 2 ] < P p 
and E[|x|| 2 ] < P s and the channel coefficients {/i™}ij,m and {^f'ji.m are i.i.d. and distributed as 
complex Gaussian CJ\f(0, 1). Each secondary receiver treats all undesired signals (interference from 
the primary users) as Gaussian interferers. Therefore, the signal-to-interference-plus-noise-ratio 
(SINR) of the n th secondary pair on the spectrum band B m is given by 

r?„E[|«| 2 ] 



SINR r 



(2) 



We define the transmission signal-to-noise ratio (SNR) by p = Throughout the paper we say 



that city and bjy are asymptotically equal, denoted by = 6jv if hniAr_ 5 . 00 = 1, and define < 



N - 

b N 

and >, accordingly. We also define the set of secondary users indices by N = {1, ... , iV}. All the 
rates in the paper are in bits/sec/Hz and log refers to the logarithm in base 2. 

2.2 Problem Statement 

Our goal is to assess the multiuser diversity gain of cognitive networks. For this purpose we identify 
M secondary transmitter-receiver pairs out of N available ones and assign one spectrum band B>m 
to each of them, such that the cognitive network throughput is maximized. We assume that all 
spectrum bands B m are of the same bandwidth. Therefore, the maximum throughput is given by 

M 



-Rmax — E 



max V log(l + SINR m , A , 
AcAf, \A\=M ^— 
m=l 



(3) 



where A m denotes the m th element of set A, for m = 1, . . . , M and the maximization is taken over 
all ordered subsets of M. 
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In order to find the optimal multiuser diversity gain in the cognitive network we first consider 
a centralized setup. We assume that there exists a central decision-making entity in the cognitive 
network, which is fully and instantaneously aware of the channel conditions of all secondary users. 
The central node solves the problem cast in ([3]) by an exhaustive search for pairing up M secondary 
users with the M available channels. For such secondary user-channel pairs we analyze how the sum- 
rate of the cognitive network scales as the number of cognitive users (N) increases. Such centralized 
setup imposes extensive information exchangqj which can be prohibitive in large network sizes. 

Next, motivated by alleviating the amount of information exchange imposed by the centralized 
setup and noting that our cognitive network is ad-hoc in nature and lacks a central-decision making 
entity, we propose a decentralized spectrum allocation scheme. In the distributed scheme each 
secondary user decides about taking over a channel solely based on its own perception of the 
network realization. We prove that the proposed distributed scheme retains the same throughput 
scaling law as in the centralized setup, i.e., is asymptotically optimal. 



3 Centralized Spectrum Allocation 

The central decision-making unit has access to all {SINR m!n } and performs an exhaustive search 
over all possible user-channel (spectrum band) combinations in order to find the one that maximizes 
the sum-rate throughput given in ([3|). In order to find the throughput scaling, we establish lower 
and upper bounds on R max , denoted by -R^ax an d -Rmax> respectively, and show that these bounds 
are asymptotically equal, i.e., R l mSkX = i?^, which in turn provide the optimal throughput scaling 
law of the cognitive network. 

We define the most favorable user of the m th spectrum band as the user with the largest SINR 
on this band, i.e., 

n* m = arg max SINR mi „. (4) 

l<n<N 

In general, it might so happen that one user is the most favorable user for two different spectrum 
bands, * / , while m ^ m! ', and as a result, these two spectrum bands cannot be allocated 

to their most favorable users simultaneously (we have assumed that each user may get access to only 
one spectrum band). Let us define T> as the event that different spectrum bands have distinct most 
favorable users i.e., no single user is the most favorable user for two distinct spectrum bands. Note 
that pairing the secondary users and the spectrum bands conditioned on the event T> is equivalent 
1 M real numbers per user. 
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to allocating each channel to its most favorable user, i.e., 

M 



E 



max log (1 + SINR m ^4 

AcAf, A\=M 



25 



E 



M 



]Tlog(l + SINR min * 



m=l 



(5) 



On the other hand, under event 25, at least one spectrum band will not be allocated its most 
favorable user and therefore we have 



E 



, ffii Sf E lo s( 1 + SINRm >^ 

AcAf, \A\=M ' 



V 



Equations (|SJ) and (jSJ) give rise to 



Ry 



E 



+ E 



< E 



< E 



M 



M 



^log(l + SINR mi<i ) 



,m=l 



max S~] log (1 + SINR m ^ 

AcJV, \A\=M ^ 
m=l 



„ M „Z>g(l + SINR fM J 
AcAA, \A =M ^— ' 
m=l 



25 



25 



P(25) 
P(X5) 



^log(l + SINR m , n;k ) 



,m=l 



-'•■max' 



Also it can be readily shown that 

Pmax > E 



„ ffii »,Zl lo g( 1 + SINR m,A m ) 
AcAf, A =M ^— ' 
m=l 



25 



P(25) 



(6) 



(7) 



(8) 



Lemma 1 P; nax and PJ^ax are asymptotically equal, i.e., R max = fi^. 

Proof: See Appendix lAl ■ 
Now, we find how PJ^ax scales as A increases. Note that the SINRs are statistically independent 
for all users and spectrum bands. The reason is that SINR m n given in ([2]) inherits its randomness 
from the randomness of g™ (fading coefficient of the channel between the n th secondary pair on 
the m th band) and {h™j}j (the fading coefficient of the channels from the j th primary user to the 
n th secondary receiver on the m th band). Since for any two different pairs of (m, n) ^ {ml , n'), the 
fading coefficients g™ and g™, refer to fading in different locations or in different spectrum bands, 



6 



therefore they are statistically independent. Similarly it can be argued that h™j and h™, j are also 
statistically independent for (m,n) 7^ (m',n'). As a result, all the random ingredients of SINR mn 
and SINR TO / jn / for (m, n) 7^ (m', re') are independent which in turn justifies the independence of the 
SINRs. Nevertheless, SINRs are not identically distributed since different users experience different 
path-losses and shadowing effects. Hence, for more mathematical tractability we build two other 
sets whose elements provide lower and upper bounds on SINR mn and are i.i.d. For this purpose we 
define 

7max = max < — > , r/ max = max 77^ 



and 7 min = min <^ — \ , r/ min = minn*. 



1 



For m = 1, . . . , M we also define the sets 5;(m) = {5;(rre, n)}^ =1 and S u (m) = {S u (m, n)}^ =1 such 
that for re = 1, . . . , N 

I „m|2 

Wmin + P s 7max Z^ jeBm |ft n j| 

and 5 u (rre,re) = — . (10) 

P»y max + P s lramZ^jel3 m \ a n,j\ 

It can be readily verified that Si(m,n) < SINR mn < S u (m,n). We use the notations <S^(m) 
and Su\m) to refer to the largest elements of the sets Si(m) and S u (m), respectively, and use 
SINRm to denote the i th largest element of {SINRm^j^Lp In the following lemma we show how 
these ordered elements are related. 

Lemma 2 For any spectrum band B m and any i = 1, . . . , N we have S®(m) < SINR$ < S$(m). 

Proof: See Appendix iBl ■ 
Now, by recalling the definition of -R r " iax given in (J7J) and noting that SINR m n ^ = SINR^ } and 
by invoking the result of Lemma [2] we get 

" M 

^log^+sfV) 1 



and iCax<E 



,m=l 
' M 



(11) 

(12) 



2 log (l +5^(771^ 
_m=l 

In order to further simplify the bounds on -R^ax gi ven m (|lip -(|12 p . in the following lemma we 
provide the cumulative density functions (CDF) of Si(m,n) and S u (m,n) ([9]) and (llOp . 
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Lemma 3 The elements of S[(m) and S u (m) are i.i.d. and their CDFs are 

Si(m,n) ~ Fi(x;m) = 1- -j^-, (13) 

p > 

"^TmaxX "I - 1 
e -x/pr) a 

and S u (m,n) ~ F u (x;m) = 1 (14) 

^"^TminX ~t~ 1 

where K m = \B m \. 



Proof: See Appendix [Cj ■ 
We denote the i th order statistics of the statistical samples Si(m) and S u (m) with parent 
distributions given in (|13p - (|14p by <S, (m) and Su\m), respectively. By denoting the CDF of 
S^' '(m) by F^ l \x; m) and that of Su\m) by Fu\x; m), for i = 1, . . . , N we have [?] 

if^m) = ^.J^fom)) (l - i^x; m))', (15) 



j=0 
i-l 



F^(x;m) = Y,( N )(F u (x;m)Y 3 (l - F u (x; m)) (16) 
By invoking the above definitions, (|lip and ()12p can be re-written as 

^max>E/ Ml + ^^rC^w), (17) 

M poo 

and i^ax < / Ml + x) dF«(*;m). (18) 
m=l J ° 

We also define 

= 1 - e~ x , (19) 

and let (x) denote the CDF of the i th order statistic of a statistical sample with N members 
and with parent distribution G(x). By using this definition we offer the following lemma which is 
a key step in finding how i?^ ax scales with increasing N. 

Lemma 4 For the distributions Fi (x;rn), F^ (x;m) and G^(x) we have 

log(l + x) dF^\x;m) 

poo 

< / log(l + p W x) dG^(x), (20) 

J 

roc 

/ log(l + x) dFi (x;m) 
Jo 



poo 

and 
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• / log(l + PWC ) dG^(x) 

'0 

K P 



log 



1 + 



P^ 



7maxP??min 



(21) 



Proof: By using the definitions of P^ (x; m) and F^ (x; to) given in ()15p -(|16 p and using the result 
of Lemma [3] we get 



F^ 1 \x;m) = ^P/(x;to) 
1 — exp 



N 



X 



PVn 



( ^7ma X a; + 1 



N 



< 



1 — exp 



G 



X 



Pllmin 

K P 



<7^7max(a;+l) 

K P 

-7max(^ + 1) 



P 



JV 



G w(. * 



max 

(x + 1) 



A' 



+ 



K P 



P 7max(z + 1) 



Now, by using (|22p and by looking at the solutions x and x' of the equations 

u = Ff\x; to), 
and it = G (1) 



Pin 



P 7max(£ + 1) 



we find that x > x', or equivalently 



G (dY (n) _^ 7max 

F«) (u;m)>± '- ir - r ^ . 

™, . T p /max 



which after some simple manipulations leads to 



log 



P ; (1) ) (u;m) + l 



> log 
- log 



r[G^) ( U ) + l 

K P 



1 + 



~pT 



7maxP^?min 



Therefore, for the lower bound on P max given in (|17p we have 

/■oo 

/ log(l + x) dPj (x; m) 
Jo 



log 



1+ F, 



(u; to) 



(in 



(22) 



9 



1 

> / log 
'0 



log 



PVmm 

K P 



(u) + l 



du 



--rra-t p 
J- ~r p TmaxP'?min 



which is the desired inequality in (j2Tj) . 

Now, note that F u (x:m) > G(—r — ) or equivalently, 

v ' v P'/max 7 



(Fi 1 )) _1 (n;m)< Wmax (GW) _1 (n). 



(23) 



Therefore, 



/■oo 
JO 



777 



log 



< / lo 
'0 



l+P??max(G (1) ) (n) 

/>oo 

/ log(l + Wmax x) dG {X) (x), 
Jo 



du 



which establishes the inequality in (|2U|) and completes the proof. ■ 
Next, by using the result of the following lemma, we establish the scaling law of -R max - 

Lemma 5 For a family of exponentially distributed random variables of size N and parent distri- 
bution G(x) ( CDF) and for any positive real number a € M+ we have 



/•oo 

/ log(l + ax)dG W (x) = log log N + log a. 
Jo 



Proof: See Appendix iDl ■ 
Now, by recalling the bounds provided in (|17|) and (| 18f) and taking into account the results of 
Lemmas [H [J] and [5] we find the optimal throughput scaling law of cognitive networks. 

Theorem 1 In a centralized cognitive network with N secondary transmitter-receiver pairs and 
M available spectrum bands, by optimal user-channel assignments, the sum-rate throughput of the 
network scales as 

R max = AfloglogiV. 

Proof: By invoking the results of Lemmas [4] and [5] on the lower and upper bounds on i? max given 
in OH) and flU) we find 



iC ax >M log log N-M log 



PVv 



K P 



7n 
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or equivalently, 



and i? max <M log log N + M log(pr] n 



lim ari fi T AT > 1 - lim 

JV->oo M log log TV N-^oo 



1. 



log 



J- "1 p TmaxP^?min 



log log N 



and 



lim 



T3U 



< 1+ lim Mfcl 
JV^oo M log log iV Vh>oo log log N 



which confirms that -R max = M log log N . This result, along with what stated in LemmaQ] concludes 
that -R max = -R max = M log log iV, which establishes the proof of the theorem. ■ 

So far, we have assumed that there exists a decision-making center that has full knowledge of 
all instantaneous channel realizations, i.e., {/i™}; and Also there is no complexity 

constraints in order to enable exhausting all the possible user-channel assignments and choosing 
the one which maximizes the sum-throughput of the network. 

The assumptions made in this section, while not being practical, are useful in shedding light on 
the sum-throughput limit of such cognitive networks. The results provided in this section can be 
exploited as the benchmark to quantify the efficiency of our distributed algorithm proposed in the 
following section. 



4 Distributed Spectrum Allocation 

In this section we offer our distributed algorithm, where each user independently of others, makes 
decision regarding taking over transmission on any specific spectrum band. We analyze the achiev- 
able sum-throughput of the cognitive network when this distributed algorithm is utilized and show 
that it is asymptotically optimal. 

4.1 Distributed Algorithm 

In order to refrain from exhaustively searching for the best user-channel matches, we consider 
assigning the available spectrum bands to only secondary users with a pre-determined minimum 
link strength. The distributed algorithm involves two major steps; normalizing the SINRs and 
comparing it with a given threshold. The underlying motivation for normalizing the SINRs is to 
balance fairness among the secondary users, in the sense that they get equal opportunities for 
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accessing the spectrum. It so happens that some transmitter-receiver pairs have a very strong link 
and some other pair a very weak link. This becomes even more likely when we have a large number 
of secondary pairs. In such scenarios if spectrum allocation is carried out merely based on the links' 
strengths, all the strong users will dominate the network and the weak users will hardly have an 
opportunity for accessing it. So for maintaining fairness, instead of comparing the links' strengths 
(or SINRs), we compare normalized SINRs. However, it should be noted that the normalization 
factors have to be designed carefully such that we do not sacrifice achieving the optimal scaling 
in favor of achieving fairness. In other words, the objective is to achieve the optimal scaling and 
fairness simultaneously. 

In the next step, the normalized SINRs are compared against a pre-determined threshold level 
and only the users with strong enough links that satisfy the threshold condition will take part in 
the competition for accessing the spectrum. 

Specifically, to each user n = 1,...,N and channel m = 1,...,M, we assign a minimum 
acceptable level of SINR, denoted by A(m, n), which is defined as follows. 

Let T(x;m,n) denote the CDF of SINR mn given in ([2]). A(m,n) is set such that we have 



Note that for any given m and n, T(x; m, n) is a non-decreasing function on [0, +00) x [0, 1] which 
ensures that there always exists a unique solution for X(m,n). Also note that as SINR min depends 
only on the incoming channels to the n th secondary receiver on the m th spectrum band, A(m, n) can 
be computed locally at the n th secondary receiver and does not impose any information exchange 
between the secondary users. It is assumed that the secondary users are aware of the number of 
secondary pairs N in the cognitive network. 

Now, each user n computes SINRi n , . . . , SINRjv/,n) normalizes them via dividing them by A(l, n) 
, . . . , A(M, n), respectively, and identifies the channel with the largest normalized SINR, and denotes 
its index by m n , i.e., 



In the next step, the n th user compares SINR^t n against X(m' n ,n) and if SINR m t n > X(m' n ,n), 
then deems itself as a candidate for accessing the channel indexed by rrin. Based on the definition in 
(|25p we define the mutually disjoint sets T-L m for m = 1, . . . , M, such that T-L m contains the indices 
of the users deemed as candidates for taking over the m th channel, i.e., 




(24) 




(25) 
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Finally, a user with its index in % m is randomly opted for utilizing the m channel. This can be 
facilitated in a distributed way via any contention based random media access method, e.g., Aloha, 
carrier sense multiple access, etc. As soon as one user takes a channel, the other users will no longer 
try to access that channel. In the following section, we analyze the sum-rate throughput scaling 
factor of the proposed algorithm. 

4.2 Sum-Rate Throughput Scaling 

We denote the sum-rate throughput by i? sum and denote the throughput of the m th channel by R m . 
Note that the construction of % m guarantees that no single user will be regarded as a candidate for 
more than one channel and also we have R sura = J2m=i ^m- By defining R m | ■ Hm as the throughput 
achieved for the m th conditioned on having users with indices in T-L m be candidates for taking over 
B m we have 



On the other hand, by noting that one member of % m will be randomly picked for accessing B 
we get 





(26) 




(27) 



From ([26]) and ([27D for any U m / we get 




(28) 



As shown in Appendix [E] we have 




in 




i=i 



(29) 



Therefore, (j28j) and (|29p together give rise to the following lower bound on R 



'hi 
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N n 



EE E 



n=l i=l 
N 

E E 

i=i 



log (l + SINR« 



p(\n m \ = n) 



n 



log(l + SINR« 



TV 

E 



P(I«J = ") 



By further denning 



and 



TV 

^-^ n \ 



t-m\ — n I ) 



AT 



logfl + SINR^ 



1=1 

we can re- write (f30j) as i? m > If we also define Qo = P{\H m \ = 0^ we get 

TV TV TV 

EQr = J p(i^i = o)+EE- p (i^ 

i=l n=i 

AT 

= E p (i^-i = n ) = 1 ' 



(30) 



(31) 



(32) 



i=0 



n=0 

which suggests that {Q™}i=o ^ s a van d probability mass function (pmf). In the sequel, we con- 
centrate on finding the scaling behavior of R l m . By using the definitions of <Sj(m) and S u (m) and 
exploiting Lemma [2j from (|32p we have 

TV 

(33) 



R l m > Q? E 

i=l 
TV 

and R l m <J2QT^ 



log (1 + 5 z W (m) 



i=l 



log fl+5«( 



(34) 



By recalling that the CDFs of <S^(m) and Su\rn) are F^(x;m) and Fu\x;m) provided in (fl~5 
and (fTUj) . respectively, (f3"3"j) and (|3"4"|) can be stated as 

TV i 

i4 > E \ + x ) dF i (t) ^ m ) 

i=l J o 
TV x 

and R l m <Y^Qf log(l + x) dF®(x; m), 
i=i ^ 

Next, for the given set of {Q^i we define 



(35) 
(36) 



TV 



Ff(x;m)±J2Q?F®(x;m), 



(37) 



i=l 
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N 



Ff(x;m)^QTF^(x;m), (38) 

i=l 
N 

and G N {x) = Y J Q?G {€) (x). (39) 



i=i 

Since {Qf 1 } is a valid pmf, F l N (x;m), F^{x;m), and G N (x) can be cast as valid CDFs. Therefore, 
([35]) and (f36|) give rise to 



I log(l + x) dF^x-m) <R l m <[ log(l + x) dF^{x;m). (40) 
Jo Jo 

The two subsequent lemmas are key in finding how i?^ scales with increasing N. 

Lemma 6 For a real variable x £ [0, 1] and integer variables N and i, < i < N — 1, i/ie function 

/(x,^£(^V-'(i-xy, 

„— n \ J / 



i=0 

zs increasing in x. 



Proof: See Appendix [Fj 

Lemma 7 For t/ie distributions F, N (x), F^(x) and G N (x) we have 

/ log(l + x)dFf(x;m) < / log(l + p7? max x) dG^(x) 
Jo Jo 



(41) 



and 



f'OO f'OO 

/ log(l + x) dff (x; m) > / log(l + P 7] min x) dG N (x) 
Jo Jo 



log 



1 + - 



(42) 



Proof: By using the definition of f(x,j) provided in Lemma [6] and recalling (|15p -(|16p we have 

F^\x;m) = f^F[(x;m),i - ij, 
F»(x;m) = /(F u (x;m),i-l), 
and G W (x) = f(G(x),i - lj . 



By following the same lines as in (|22j) we also have 



(^ + ^7max(x + l) 
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and consequently by applying Lemma [6] and using the definition in (|37p - (|39p we have 
Ff(x;m) 

N N 

= Y, QT F ® (x;m) = J2QTf(F l (x;rn),i- l) 

i=l i=l 

F'/min -<s /I 



N / 

i=i v 

= E^G«(-f- + ^7 m a X (x + l)) 

. , Vf'/mm r s / 



= G I 1 — — 7 max (x + 1) 

By following a similar approach as in Lemma 01 and f|22|) - ()23|) the inequality in (|42|) can be estab- 
lished. Proof of (|41|) follows a similar line of argument. ■ 



Lemma 8 For a family of exponentially distributed random variables of size N and parent distri- 
bution G(x) (CDF) and for any set of {Q^^fLi such that YliLoQT 1 = 1j if the condition 

Jim ^=i ^ =0 (43) 
is satisfied, then for any positive real number a £ M+ we /iaue 

/>oo 

/ log(l + ax)ciG Ar (x) = log log TV + log a. 
Jo 

Proof: See Appendix [Gj ■ 
By using the results of the Lemmas [7] and [8] we offer the main result of the distributed algorithm 
in the following theorem 

Theorem 2 The sum-rate throughput of the cognitive network by exploiting the proposed distributed 
algorithm scales as 

Rsum = M log log N. 



Proof: We start by demonstrating that the set {Q™} as defined in (|3ip fulfils the condition ([33 
of Lemma From (|3ip we have 

N N N 



i=l i=l n=i 

N 



j=l n=i 

AT 



n 

n=l i=l 
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N , 1 
—P{\n m \ = n 

71=1 
AT 

n 



71=1 
1 ^ 

+ 9E P (l^™l = n ) - o P (l^-l = °)- ( 44 ) 



n=0 



= 1 

Note that |% m | has compound binomial distribution with parameters {p(m,n)}^ =1 |10j . where 
p(m, n) denotes the probability that the m channel is allocated to the n th user. Therefore, 
according to the properties of compound binomial distributions we have jlO] 

N N 

^nP(\H m \=n) =E[\'H m \]=J2p(m,n). (45) 

n=l n=l 

From (El) and (USD we get 



JV / N 

7=1 \n=l 



p(m, n) + 1 . 



On the other hand, the probability that any specific user n can be a candidate for taking over any 
of the M channels is 

, x a / f SINR m , n ) 

w(n) = max < — ^ > > 1 

\ m [ A(m, n) J 

M 

= 1- J] P ( SINR ™,n < X (m,n) 

771=1 

M ( \ ( 1 \ M 

= 1- [[ T[X(m,n);m,nJ = 1- ( 1- — J . 

777=1 " v ' ^ ' 

= 1-1/AT 

Therefore, since X^m=i P( m > n ) = w ( n )> f° r an w-,^ we have p(m,n) < co(n). Hence, 



N 1 

^iQr<-(iVa;(„) + l) 



7=1 

On the other hand, 

lim Nuj(n) = lim v = N! = M. 

Therefore, YliLi i-QT — \iM 1) an d the set {Q™} satisfies the condition in Lemma IH1 Therefore, 
Lemmas [7] and [8] together establish the following 



poo 

/ log(l + x) dFj v (x;m)<loglogiV + log(pr? max ), 
J o 
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Figure 1: Sum-rate throughput versus the number of secondary users for M = 1, . . . , 4, p = 10 dB, 
and the number of the primary users K m = 4. 



POD 

and / log(l + x) dF t N (x; m)> log log N 
Jo 

-log 



K P 

P 
1 s 



The two inequalities above, in conjunction with PU|) and noting that -R sum = Ylm=i provide 

1 



Mlog log N-M log 



+ 



K P 

J1 m J p 



7n 



Rsum 

< M log log iV + M log (pr/n 



which concludes the desired result. 



4.3 Simulation Results 

The simulation results in Fig. [1] demonstrates the sum-rate throughput achieved under the cen- 
tralized setup given in (j3J) and the distributed setup given in Section 14.11 We consider a primary 
network consisting of 4 users and look at the throughput scaling for the cases that there exist 
M = 1, ... ,4 available spectrum bands to be utilized by the secondary users. We set all path-loss 
terms {r/j} and {"fij} equal to 1 and find the sum-rate throughput as the number of secondary 
users increases. As shown in Fig. [lj as the number of secondary users increases, the sum-rate 
throughput achieved by the centralized and distributed schemes exhibit the same scaling factor. 



18 



1 — ' — ■ — 1 
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Number of Secondary Users (N) 

Figure 2: A(m, n) versus the number of secondary users for M = 4 available spectrum bands and 
K m = 4 primary users. 

Note that what Theorems 1 and 2 convey is that the ratio of i? max and -R S um in the centralized and 
distributed setups, respectively, approaches to 1 as N — > oo, i.e., 



iV->oo i? sum 

which does not necessarily mean that i? max arid i? S um have to coincide. As a matter of fact, as 
observed in Fig. 1, there is a gap between i? max and -R SU m> that according to the results of Theorems 
1 and 2 must be diminishing with respect to -R max and -R S um such that we obtain the asymptotic 
equality i? max = R S um = M log log N. This gap accounts for the cost incurred for enabling 
distributed processing in the distributed spectrum access algorithm. 

The throughput achieved under the distributed setup uniformly is less than that of the central- 
ized setup. This is justified by recalling that the centralized scheme finds the best secondary user 
for each available spectrum band, whereas the distributed network finds all the secondary users 
whose quality of communication on a specific channel satisfies a constraint (A(m, n)) and among all 
such secondary user one is randomly selected to access the spectrum band. This, not necessarily 
guarantees finding the best user for each available spectrum band and as a result leads to some 
degradation in the sum-rate throughput. 

Finding the metric \(m, n) as defined in (|24p is the heart of the distributed spectrum allocation 
algorithm. As it is not mathematically tractable to formulate the CDF of SINR m)n (note that 
F[(x; m) and F u (x; m) are only the CDFs of the lower and upper bounds on SINR m>n ), we are not 
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able to find a closed form expression for A(m,ra). However, by solving (|24p numerically, we provide 
the following two figures which are helpful in shedding light on how A(m, n) varies with other 
network parameter, i.e., primary and cognitive network sizes as well as the number of available 
spectrum bands. 

Figure [2] demonstrates the dependence of A(m, n) on the transmission SNR denoted by p. It is 
seen that X(m, n) monotonically increases with p. Intuitively, as p increases, the users are expected 
to have more reliable communication and as a result the algorithm will impose more stringent 
conditions on the secondary users for considering themselves as a candidate for accessing any 
specific spectrum band. More stringent conditions will translate to having higher values of A(m, n) 
such that the condition in f)24f) is satisfied. 

The numerical evaluations provided in Fig. [3] show that A(m, n) increases as the size of the 
primary network decreases. Again as in Fig. [21 this is justified by noting that smaller number 
of primary users leads to less interference from the primary network to the cognitive network and 
thereof, more reliable secondary links. Thus, decreasing the primary network size again requires 
more stringent conditions to be satisfied for a secondary user to be deemed as a candidate for 
taking over a spectrum band, which in turn results in an increase in X(m, n). It is noteworthy that 
the choices of the thresholds given in (I24D have been heuristic choices that satisfy all the desired 
properties (optimal scaling as well as fairness and limited information exchange as discussed in 
Section [5]). Nevertheless, we cannot prove that these are the only choices of the thresholds and it 
might be possible to find some other threshold settings that satisfy all these conditions and yet do 
not depend on the size of the primary network. Hence, while the scaling of the sum-rate throughput 
does not depend on the size of the primary network, the choices of the thresholds for achieving this 
scaling in the distributed algorithm do depend on the size of the primary network. 
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Figure 3: A(m, n) versus the number of secondary users for different sizes for the primary network 
K m = 1,...,4 and p = 10 dB. 

5 Information Exchange and Fairness 
5.1 Information Exchange 

In the distributed algorithm we assume that the n th secondary receiver measures {SINR mjn }^ =1 cor- 
responding to different spectrum bands, selects the largest one, whose index is denoted by m^, and 
compares it against a pre-determined quality metric A(m, n). If SINR m t n > A(m, n), this secondary 
receiver should notify its designated secondary transmitter to participate in a contention-based com- 
petition for taking over channel rati- Such notification requires transmitting logM information bits 
from the secondary receiver to its respective secondary transmitter. 

Although not all of the secondary pairs will be involved in such information exchange, it is 
imperative to analyze the aggregate amount of such information for large networks (N — > oo). In 
the following theorem we demonstrate that for the choice of A(m, n) provided in (I24D . the asymptotic 
average amount of information exchange is a constant independent of iV and therefore does not 
harm the sum-rate throughput of the cognitive network. 

Theorem 3 In the cognitive network with distributed spectrum access, when X(m,n) satisfies 

T^X(m,n);m,nj = 1 - — , 

the average aggregate amount of information exchange between secondary transmitter-receiver pairs 
is asymptotically equal to MlogM. 



21 



Proof: 
1-(1 
Rie, is 



As stated earlier, the probability that the a user satisfies the \(m,n) constraint is co(n) = 
- l/N) M . Therefore, the average aggregate amount of information exchange, denoted by 



R ie = lim N Lo(n) log M 



iV-j-oo 




1 \M-1 



M log M, 



(46) 



(47) 



where for the transition from (|46p to (|47p we have used the L'Hopital's rule. 



5.2 Fairness 

In general, opportunistic user selections might lead to the situation that the network be dominated 
by the secondary pairs with their receiver far from the primary users so that see less amount of 
interference from them, or by those pairs where the transmitter and the receiver are closely located 
and enjoy a good communication channel. 

In despite of these facts, we show that in our network, by appropriately choosing {A min } we can 
provide equiprobable opportunity for all users to access the available spectrum bands. This can 
be made possible by enforcing more stringent conditions (higher \(m,n)) for the users benefitting 
from smaller path- loss and shadowing effects. In the following theorem we show that by the choice 
of A(m,n) provided in (|24p all the users have the equal opportunities for accessing a channel. 

Theorem 4 In the cognitive network with distributed spectrum allocation, when A(m,n) satisfies 



all the users have the same probability for being a allocated a channel. 

Proof: As shown earlier, the probability that user n satisfies the SINR constraint A(m, n) is 



6 Discussions 

6.1 Impact on the Primary Network 

In cognitive networks with underlaid spectrum access, the secondary and primary users may coexist 
simultaneously. Therefore, in order to protect the primary users, the secondary users must adjust 



T[\{m,n)\m,n) =1- — 



uj{n) = 1 — (1 — 1/N) M , which is the same for all users. 
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their transmission power such that they operate within the tolerable noise level of the primary 
users and thereof do not harm the communication of the primary users. Hence, it is imperative 
to investigate whether such power adjustments affect the achievable throughput scaling in the 
centralized and distributed setups. 

According to theorems [T] and [2 the sum-rate throughput of the cognitive network scales as 
M log log iV which does not depend on the SNR or the transmission power of the secondary users. 
Therefore, irrespective of the transmission policy and power control mechanism (i.e., for any ar- 
bitrary SNR or transmission power), the secondary users achieve the scaling law of M log log N. 
Hence, deploying any power management mechanism of interest along with the proposed spectrum 
access algorithms, does not harm the optimal scaling. 

6.2 Distributed Algorithm 

For implementing the distributed spectrum access protocol there are two major steps involved. 
First the random selection of a user out of the candidates for taking over a specific spectrum band. 
For randomly selecting a user out of the set of users in % m to access the m th spectrum band, one 
distributed approach is to equip all the users with a backoff timers. Then when a user learns that 
it is a candidate for accessing the m th spectrum band with run the backoff timer with an initial 
random value. The first cognitive pair whose backoff timer goes off will take over the channel and 
with a beacon message can notify it to the rest of the network. 

Secondly, the distributed algorithm requires some secondary receivers transmit log M informa- 
tion bits to their respective transmitters. Transmitting log M information bits requires a very low 
rate communication. An appropriate approach for such communication rate is to deploy ultra-wide 
band (UWB) communication between a secondary transmitter and receiver pair. This will allow 
the secondary users to communicate the low-rate information bits well below the noise level of the 
primary users. It is noteworthy that the cognitive radios are often assumed to be equipped with 
wideband filters which enable them to transmit and receive in a wide range of frequency spec- 
trum. This feature of the cognitive radios provides an appropriate context for implementing UWB 
communication. 

6.3 Networks of Limited Size 

Practical networks do not have large enough number of users to fully capture the multiuser di- 
versity gain (double- logarithmic growth of capacity with the number of users). Therefore, due to 
such degradation in multiuser diversity gain, the network cannot support the throughput expected 
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in theory. Therefore, for practical networks we have only upper bounds on the actual sum-rate 
throughputs. Knowing such upper bounds help to find in an insight about what to expect from 
the cognitive networks at the design state. 

Although the result hold analytically for only N — > oo, from the simulation results in Fig. 1 we 
observe that as low as N = 50 (which is around the point that we start to observe steady increase 
in the throughput) secondary users are enough to start observing the multiuser diversity gain. This 
is not far from the size of practical networks. 

7 Conclusions 

In this paper we investigated the multiuser diversity gain in cognitive networks. We first obtained 
the optimal gain achieved in a network with a central authority and show that the gain achieved 
in such cognitive networks is similar to that of the interference-free networks, i.e., the network 
throughput scales double-logarithmically with the number of users. Then we proposed a distributed 
spectrum access scheme which is proven to achieve the optimal throughput scaling factor. This 
scheme imposes the exchange of log M information bits per transmitter-receiver cognitive pair for 
some pairs, and no information exchange for the others. The other specification of the distributed 
algorithm are that the network-wide average aggregate amount of information bits it requires is 
asymptotically equal to MlogM, and it ensures fairness among the secondary users. 

A Proof of Lemma [I] 

We equivalently show that limjv^oo P(T>) = 1. An intuitive justification is that if we put the SINRs 
in an M x N array, and locate the maximum element of each row, event T> occurs when no two 
such maximum are located in the same column. Therefore, as the number of the columns increases, 
in the asymptote of very large values of N, the even T> must occur with probability 1. 

For set C C M = {1, . . . ,M} such that \C\ > 2 we define P(C) as the probability that the 
spectrum bands with indices in C have the same most favorable user. Therefore, we get 



P(C) 



A 



P { V m, m' G C, n* 



n m'} 



N 



P {V m G C, n, 



n n* = n 



} 



71=1 



xP{n* m = n) 



N 




n=l mgC 
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xP(n* m =n) (48) 

N 

= e n n p ( sinr -,- > siNR mV ) 

n=l mgC n^n' s v ' 

= q(m,n,n') 

xP(n* m = n) 

< N(q m ^ N ~ l \ (49) 

where g ma x = iiiax m n n / g(m, ra, n') and (}4"8|) and ([4*9]) hold due to the statistical independence of 
the elements in {SINR m n }, and . Therefore we have 

p(v) = 1- Yl p (°) 

CCM,\C\>2 
M 

= i-E E ^) 

m=2 CCM, |C|=m 
m=2 s 



m(7V-l) 

max J 



as jv->oo 

= 1, 

which completes the proof. Note that q max is a function of {r^} and {Ajj} and does not depend on 
N. 



B Proof of Lemma [2] 

We first show that for any i = 1, . . . , N , 5^(m) < SINR^ • For i = 1 we have 

SINR.W = maxSINR mn > max 5/ (m,n) = 5, {1) (m). 

n n 

For any i = 2,...,N, from the definition of SINRm and Sm it can be deduced that each of 
the (JV - i + 1) terms SINR$, . . . , SINR^ } is greater than one corresponding element in the set 
Si{m). Therefore, there cannot be more than (i — 1) elements in Si{m) which are all greater than 

sinr#,...,sinf$°. 

Now, if 5^(m) > SINRm , then all the i terms 5^(m), 5^(m), . . . ,5^(m) should be greater 
than all the (N-i + 1) terms SINR^, . . . , SINR^ } . Therefore, we have found i elements in Si{m) 
that are all greater than SINRm , ■ ■ ■ , SINR^ and this contradicts with what we found earlier. 
Hence, we should have sf\m) < SINR$. 

By following the same lines, we can show that also for i = 1, . . . , N, we always have SINR^ < 
Su (jn) , which concludes the proof of the lemma. 
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C Proof of Lemma [3] 

Let Y = \g™\ 2 , which has exponential distribution with unit variance. Also define Z = J2j=i \ h n,j\ 2 
which is the the summation of K m independent exponentially distributed random variables each 
with unit variance, and thereof has a Gamma(if m , 1) distribution. By denoting the probability 
density functions (PDF) of Z and Y by 



My) 

and f z (z) 



*gK>m 1 g z 



{K m -l)V 

the PDF of Si(m,n) = i/ p??min+i yp 37maxZ , denoted by f s (x) is 

/•oo 

fs(x) = I fs\ z(x | z)f z (z)dz 

1 I Ppln 







V PVmin Ps 



fp7m ° zi| 



g Km 1 g Z 



e 



-x/pVmin 



max 



— ;,t p 

PVm'm *s 

where the last step holds as e~ u u M = Ml. Therefore The CDF is 



F[(x; m) = 1 



s K ■ 

V 



F u (x;m) can be found by following the same lines. 



D Proof of Lemma [5] 



We start by citing the following theorem. 



Theorem 5 jilt Theorem 4] Let {X n }„ =1 be a family of positive random variables with finite 



mean /ijy and variance a 2 N , also /ijy 



oo and 

(J-N 



as N — > oo. Then, for all a > we have 



E 



log(l + aX N ) = log ( 1 + aE[Xjv] 



2G 



Consider the set of random variables {Yi, . . . , Y/v} where Yj ~ G(y) and define X^ = y^ n ~ 1+1 \ 
where is the order statistic of the set {Yi, . . . , Yjv}; hence, ATtv ~ G^(x). Therefore, as 
provided in [?, Sec. 4.6] 

N 



UN 



^■2 A 



e[x n] = y: \ 

n=l 

E[\X N - m \ 2 ] = ^2^, 

71 
n=l 

which confirms that for finite N, fJ,N is also finite. Also as shown in [3], fi^ = logiV and af^ = 
from which it is concluded that — — > as N — > oo. Therefore the conditions of the above theorem 
are satisfied and we have 

roo 

log(l + ax)G (1) dx = E log(l + aXjv) 



= log (l + aE[X N ] 

= log(l + a log N) 
= log log N + log a 



which is the desired result. 



E 

By using the following lemma, we further find a lower bound on R m which will be more mathe- 
matically tractable. 

Lemma 9 For a continuous random variable X, increasing function g(-) and real values b > a 



E 



g(X) \ X>b > E g(X) \ X > a 



Proof: See Appendix IHl 

(i) 

By recalling the definition of SINRm we have 



E log(l + SINR r 
Y E log(l + SINR r 



SINR mji > X(m,i) 



w , SINR m , 4 SINR m /,i 
;Vm m : — > - 

\{m,i) A{m,i) 



(50) 
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>E E 
>E E 

i€H m 

= E E 



log(^l + SINR r 
log(l + SINR r 
log(l + SINR r 



SINR m ,i > X(m,i) 

SINR m j > min A(m, j) 
j 

SINR mj j > minA(m,j) 

j 

;Vl(£H m : S\NR hn < minA(m,j) 



(51) 
(52) 



\Hm\ p 

]TE log(l + SINR# 

|«m| 

'log(l + SINR# 



SINR^ > minA(m,: 



(53) 
(54) 

(55) 



>E E 

j=i 

where (|50p is obtained by replacing "H m by an equivalent representation. Transition from (|50p to 
(|5ip holds by applying Lemma [9] for 6 = SINR m / j • j^^k and a = for all to' ^ to. Transition to 
(|52p is again justified by using Lemma[9j Due to the statistical independence of SINR m j and SINR m / 
for to G % m and / ^ "H m the additional constraints imposed in (|53p do not result in any changes. 
The conditions in (|53p are equivalent to having the |"H m | largest components of {SINR^n}^^ be 
greater than mini<j<Ar A(i), which is mathematically stated in (f54|) . Finally, (J55j) holds by applying 
Lemma [9] one more time. 



F Proof of Lemma [6] 

By the expansion of ( x + (1 — x) J we have 



iV 



j=i+l VJ 



v 



A/ 



j=i+i 

N-(i+l) 



k 



fc=0 

l-/(l-x,iV-i-l) 



A/ 



1 - x) JV - fe x fe 



where it can be concluded that f r (u,i)\ _ = f'(u,N — i — 1)| _ 1 _ • So it is sufficient to show 
that f'(x, i) > for x < \ and for all i = 1, . . . , N — 1. For this purpose we consider two cases of 
*< Lfj andi> Lfj. 
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Case 1: i < |_f _ 



f'(x,i) 



^-i-\\-x) 



3=0 

N 



jx N ~ j (l-x) 



i=o 



N 



x N -i-\\- x )i- 1 N(l-x)-j 



where since < j < i it can be shown that for x < | 



2V(1 -z) - j > iV(l -x)-i 

N 

>N(l-x)-- 



N 
~ ~2 
> 0. 



(1 - 2x) 



(56) 



Case 2: i > [f J 

Define ctj = 1 — ^<5(|_irJ ~~ i)' where <5(-) is the Dirac delta function. Therefore, we get 



f(x,i) = f(x,N-i-l) 



LfJ 



N 



+ S a i x^Cl - xy + x j (l - x) N ~ j 



For x < \ we get 



/'(x,i) = /'(x,JV-i-l) 



LfJ 



TV 



J=iV-i v J 



N - j - Nx 



+ x j ~ l {l - x) 



N-j-l 



j — Nx 



>x N -i- 1 {l-x)i- 1 
> f'(x,N-i-l) 



<LfJ 




> 0. 



(57) 
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From (|56p and (|57p it is concluded that for x < ^, f(x,i) is an increasing function of x, which 
completes the proof. 

G Proof of Lemma [8] 

This proof follows the same spirit as the analysis provided in [JJ. However, due to some differences 
in our setting, we provided an independent treatment. 

For any given number of users N, we define a random variable Xn, distributed as X/v ~ G N (x) 
and also for j = 1, . . . , N we define 



A 

m = 



oo 

X I 



dG {j) { 

and <Ty) = y dG^\x), 

/•oo 

and /xjv = E[Xtv] = x dG N \x) 

Jo 

AT oo JV 
„■ i „■ i 



3=1 ^ 3=1 

As given in [?, Sec. 4.6] and discussed in details in [3J, for ordered exponentially distributed random 
variables we have 

2 



a N < 2 + 2/i (1) ( - n N J , (58) 
1 

l(N + 



and logAr + C+ ^ </i ( i) <logiV + C + ^, (59) 



and therefore, 

M(l) = log iV, 

where £ ~ 0.577 is the Euler-Mascheroni constant. Also 

N 



M(i) - lo g ( X] ^ ) - C - 0.5 < /iAr < 



' 3= 

By taking into account the constraint in (03]), as N — > oo 



log^EjLiJ'Qi) -C-o.5 



log iV (60) 
Equations (f5T?j) and (f6"Uj) together show that 

M(i) = /^JV = log X, (61) 
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which also implies that /zjv —too. Taking into account (|58p and (|6ip we also conclude that 
liniAr^oo Z£L = and therefore the conditions of Theorem [5] are met. Hence, from Theorem [5] 



log(l + ax) dG N {x) 



= E log(l + aXjv) 

= log (l + aE[X N ] 

= log (l + a/j, N 

= log log iV + log(a) . 



H Proof of Lemma [9] 



E 



g{X) \ X>b 



g{x)f x \x>b{x) dx 



1 



> 



Pr(X > 6) Jb 
Pr(X > b) 
~ Pr(X > a) 



1 



+ 



> 



+ 



Pr(X > a) Jb 

9(b) 
Pr(X > a) 
1 

Pr(X > a) 
1 

Pr(X > a) 
1 



9{x)fx(x) dx 

9(b) 
g(x)fx(x)dx 
Pr(a < X < b) 
g(x)fx(x)dx 

g(x)fx(x) dx 
g(x)fx(x)dx 



g(x)fx\X>a(x) dx 

E g(X) | X > a . 
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